Longitudinal muon spin relaxation in high purity aluminum and silver 
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The time dependence of muon spin relaxation has been measured in high purity aluminum and 
silver samples in a longitudinal 2 T magnetic field at room temperature, using time-differential 
/i + SR. For times greater than 10 ns, the shape fits well to a single exponential with relaxation rates 
of Aai = 1.3 ± 0.2 (stat.) ± 0.3 (syst.) ms" 1 and A Ag = 1.0 ± 0.2 (stat.) ± 0.2 (syst.) ms" 1 . 



PACS numbers: 76.75.+i 



I. INTRODUCTION 



When positive muons are stopped in a high purity 
metal with the muon spin along the direction of an ap- 
plied magnetic field, the muons can be depolarized by 
interactions with nuclear dipole moments, conduction 
electrons, and paramagnetic impurities. The form of 
the resulting depolarization has been studied using the 
/t + SR technique,— but has not been quantified to the level 
needed for the TWIST experiment, 2 in which the abso- 
lute polarization of the muons must be known with high 
precision. TWIST used stopping targets of aluminum 
and silver with purity greater than 99.999% immersed 
in an external 2.0 T longitudinal magnetic field at room 
temperature. Because the TWIST drift chamber detec- 
tor system recorded ionization from incoming muons as 
well as positrons with a relatively slow drift gas, the anal- 
ysis discarded decay positrons within 1.05 /is of the muon 
stopping time due to possible time overlap of muon and 
positron ionization. Therefore, a time differential /i + SR 
experiment was undertaken to make a precise measure- 
ment of the relaxation rate, especially for decay times 
below 1.05 /is. 



II. DEPOLARIZATION MECHANISMS FOR 
STOPPED MUONS 

After motional thermalization, positive muons are lim- 
ited to interstitial and substitutional sites. When nearly 
thermalized, a muon in a good metal attracts a screening 
charge of conduction electrons^ At room temperature, a 
muon is usually diffusing between energetically allowed 
sites before decaying. The conduction electrons in alu- 
minum and silver efficiently screen the ionic potentials 
over all distances greater than the Fermi- Thomas screen- 
ing length. 

Muons can become trapped at crystal defects, of which 
there are a wide variety; the most common originate 
from the manufacturing process, such as when a metal 
is cold-rolled to produce a thin foil. Trapping by such 



defects makes the muon mobility strongly sample depen- 
dent. The defects are enhanced by quenching, and can 
be diminished by annealing. 

In slowing down from high energy, the muon itself can 
cause lattice defects as it imparts recoil energy to lat- 
tice ions in ~ 10 _17 s and the lattice distributes this en- 
ergy to neighbouring atoms in ~ 10 _12 s. 4 A nucleus can 
be knocked out of its lattice position into an intersti- 
tial site, leaving a vacancy (Frenkel defect). However, 
these vacancies are unlikely to affect the muon's diffu- 
sion, since the muon thermalizes ~ 1 /zm from the last 
defect introduced^ 

The magnetic field experienced by a stationary muon 
due to the magnetic dipole moments of nuclei and lattice 
impurities can be modelled as static, isotropic, and Gaus- 
sian. If the muon samples a new field taken at random 
from the same distribution every time it "hops" to a new 
lattice site, the resulting depolarization is given by£ 



Pu(t) = P^(0)exp 



2A 2 



[exp(-ut) - 1 + vt] \, (1) 



and is valid for i//A sufficiently large, where A is a mea- 
sure of the magnetic field distribution and 1/v is the 
mean time between muon hops. Each field component 
is presumed to have an independent Gaussian distribu- 
tion P(Biocai) ~ exp [-B 2 ocal /(2A 2 / 7 2 )] , where 7(U is the 
muon's gyromagnetic ratio and A/7^ is the standard de- 
viation of a Gaussian distribution of magnetic fields. 

If an external field B cxt is now applied in a direc- 
tion transverse to the muon polarization, the muon spins 
precess and are depolarized according to the Abragam 
formulaf£ 



cos (tu^i) 



P„(t) = P M (0)exp j-^L [exp(-*/t) - 1 + irf] J 

(2) 

where uj^ = 7^-Bext- In the "motional narrowing" limit, 
the muons move quickly so that v is large, exp (—ft) — > 
and the envelope of Eq. ((2J approaches an exponential 
time dependence. In the static limit, the envelope instead 
approaches a Gaussian time dependence. 
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If a longitudinal field B$ is applied, the static relax- 
ation function becomes 8 - 

2A 2 

P»(t) = 1 2" [ X ~ ex P (~3 A2 * 2 ) cosw i] 

2A 4 /"* / 1 \ 

H 5- / exp --A 2 r 2 sin uj t dr , (3) 

wo io V 2 J 

where uio = 7 M -Bo ~ 272 MHz, and the longitudinal 
field is seen to suppress the depolarization due to nu- 
clear dipole moments. The largest observed nuclear dipo- 
lar field on a muon in a crystal is = 0.47 mTp 
while the applied longitudinal field in this experiment 
was Bo = 2 T, so that (2A 2 /w 2 ) < 1(T 7 . Of course, in 
our experiment the muons are not static. The cumlativc 
loss of polarization increases with increasing v as long 
as v <C ujq- For v 3> luq "motional averaging" sets in 
and the relaxation rate goes down again; this is known 
as the "Ti minimum" effect^ So the fastest relaxation 
would be for v ~ 100 /is^ 1 ; assuming muons lose up to 
10~ 7 of their polarization at each hop, in that case the 
relaxation rate is < 0.01 ins -1 . At room temperature 
in aluminum 1 ^ v ks 10 11 s _1 , which is far past the "Tl 
minimum" . Depolarization by nuclear dipole moments is 
therefore negligible for this experiment. 

Korringa relaxation— caused by a hyperfine contact in- 
teraction between the muon spin and fluctuating conduc- 
tion electron spins is the most likely cause of any muon 
depolarization in this experiment. The net hyperfine cou- 
pling is an average over all electron spin orientations . 13 i 14 
The resulting exponential relaxation rate is given by 

P M (i) = P M (0)exp(-Ai). (4) 

where A is proportional to K^T. The muon Knight shift 
= (uj — ujq)/luq is the fractional muon frequency shift 
due to the same hyperfine interaction. Korringa relax- 
ation has been observed in several non-magnetic metals 
(lead, cadmium, zinc, copper), where the muon relax- 
ation rates increase linearly with temperature and are 
robust to field changes in the range 0.01 T to 0.20 T^ 
as expected^ 

As we have seen, relaxation of muons by nuclear dipole 
moments is heavily suppressed by a longitudinal mag- 
netic field. The electronic dipole moments of param- 
agnetic ions are much larger and can produce fields up 
to 0.1 T at a distance of one lattice spacing j 16 ' 17 Thus 
they are much more difficult to decouple from the muon 
and may cause exponential relaxation of muons in a 
crystal with significant concentrations of paramagnetic 
impurities il i 16 ] 18 

In summary, muon spin relaxation due to nuclear 
dipole moments exhibits an exponential form when the 
muons hop rapidly, as expected at room temperature. 
Non-exponential static dipolar relaxation may occur if 
the muons become trapped at defects or vacancies, which 
is possible in our aluminum foil since it was not annealed, 
but any such static relaxation due to nuclear dipole mo- 
ments is heavily suppressed by the presence of a strong 



longitudinal magnetic field. In the silver foil, trapping 
at paramagnetic impurities cannot be excluded. The ap- 
propriate form is still exponential, as long as the muons 
diffuse sufficiently fast to find the impurities promptly. 
This may be the case in an annealed sample such as our 
silver foil. For both foils, Korringa relaxation is expected 
to be the dominant mechanism for depolarization. 



III. PREVIOUS MEASUREMENTS 

We were unable to find any measurements from /i + SR 
taken under conditions comparable to those of TWIST. 
Depolarization in aluminum and silver from nuclear 
dipole moments has been measured in /i + SR experi- 
ments, using a transverse magnetic field arrangement. 
There are more studies on aluminum since its nuclear 
dipole moment is about 35 times larger than silver. 
Even with its large dipole moment, high purity alu- 
minum leads to almost negligible depolarization down to 
1 K»i2r— Consequently, a measurable depolarization is 
typically seen only in samples that have been doped with 
impurities 4 i 24 ' 25 

There is a contradictory measurement^ which used 
aluminum and silver targets of 99.99% purity in a trans- 
verse field arrangement, at room temperature, and ob- 
served a Gaussian form for the depolarization in alu- 
minum. They explained this anomalous result as due 
to muons trapping in defects, which originated from the 
cold-rolling during manufacture of the foili^ 7 - 

The result from an intermediate phase of TWIST ^ 
which used a 2 T longitudinal magnetic field, is Aai = 
1.6 ± 0.3 ms _1 . One other measurement^ 8 - was available 
prior to this experiment, where a high purity aluminum 
sample was measured in a 1.1 T longitudinal field. The 
result, based on four points in the time range of 0.1 to 
10 /is, is Aai = 0.43 ± 0.34 ms _1 . It is difficult to com- 
pare the results from these two experiments due to dif- 
ferences in magnetic fields, limited statistical precision, 
and suspected systematic uncertainties in determining 
A. Comparing results between different experiments has 
the additional complication that if the depolarization is 
dominated by small impurities and/or defects introduced 
during manufacturing, then the result is expected to be 
dependent on the sheet from which the sample is taken. 



IV. EXPERIMENTAL DETAILS 

In December 2006 we acquired ^ + SR data using the fa- 
cilities of the TRIUMF Centre for Molecular k Materials 
Science. The M20B beamline, which was a dedicated sur- 
face muon channel, 29 included collimators for the muons 
and a DC separator that removed positrons by velocity 
selection. The separator also rotated the muon spin by 
11°, creating a small transverse polarization component 
and reducing the longitudinal polarization by 2%. The 
muon rate was typically 20-40 kHz. The muon target 
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and positron counters are shown schematically in Fig. [T] 
Muons could stop in the sample under investigation, one 
of the scintillators, or an Ag mask of purity 99.99%. 
Most data were taken with a TM2 scintillator of nominal 
thickness 254 /im, but data were also taken with nominal 
thicknesses of 127 and 508 /im. The backward (B) and 
forward (F) positron counters were approximately 0.6 cm 
thick. The backward counter was a disc of about 8.0 cm 
diameter, with a 2.5 cm hole for the muon beam. The 
F counter was tubular, with inner diameter of 10.5 cm 
and length 35.5 cm. The target module was placed in the 
HELIOS superconducting solenoid that supplied a longi- 
tudinal 2 T magnetic field over the stopping target. The 
target was at room temperature, and the region around 
the target was evacuated. 

The aluminum stopping target was purchased from 
Goodfellow Corporation,— who specified the typical im- 
purities as 0.3 ppm of Cu, 0.3 ppm of Fe, 1.2 ppm of 
Mg, and 0.8 ppm of Si. In aluminum, the muons are 
not trapped by impurities above wlOO KM High pu- 
rity aluminum has been studied under annealing and 
quenching^ over a temperature range of 19 to 900 K: 
most defects were found to be absent after allowing the 
quenched sample to reach room temperature. The sil- 
ver stopping target was purchased from ESPI Metals^ 
who specified the typical impurities as 2 ppm of Fe, < 2 
ppm of Bi, 0.6 ppm of Cu, and 0.6 ppm of Pd. It was 
annealed in an inert argon atmosphere, after machining. 
In silver, there is evidence that room temperature trap- 
ping at impurities can occur ji£ but that would not be of 
consequence for our annealed sample. 

The data acquisition system recorded the positron time 
of arrival relative to that of the muon, in two pairs of his- 
tograms, each with 19200 channels of width 0.78125 ns. 
Data were recorded for the B and F counters, in par- 
allel for muons stopping in the Ag mask and for those 
stopping in the target, separated by the absence or pres- 
ence of a signal in TM2 (see Fig. [1} . There was a back- 
ground of muon decays in or near the collimator, scat- 
tered positrons, and muon decays between the muon de- 
tector and the target. The electronics removed the back- 
ground from muon pile-up, where a second muon enters 
the sample region during the same 15 /is long data gate as 
the original muon. The positron gate was open for 1 /is 
before the muon trigger, so that events were recorded 
that could not possibly be caused by decay of the de- 
tected muon. The aim was to fix the background to the 
average of the number of decays recorded in this 1 /is 
interval. 

The stopping targets are listed in Table HI For sets E- 
K, the target consisted of two layers, one of metal from 
the same foils as the TWIST targets, the other a cali- 
bration sample (CS). The orientation of the layers could 
be reversed to expose either one to the stopping muon 
beam. The metal portion was either 10 folded layers of 
Al or 6 folded layers of Ag, both of which were sufficiently 
thick to stop all the muons in the surface muon beam. 
The layers were bound together with 3.6 /xm thickness 



NOT TO SCALE 




vacuum window 



B = backward counter, disc 
F = forward counter, cylinder 
TM1,TM2 = muon scintillators 

FIG. 1. Schematic for the experiment. A collimated muon 
beam enters from the left and muon decays are recorded in 
one set of histograms if they are detected only by the first 
scintillator (TM1), stopping mostly in the Ag mask, or in 
another set of histograms if they are recorded by both scintil- 
lators (TM1 and TM2), mostly stopping in the target. The 
backward positron counter (B) is protected from the muon 
beam by a collimator. Materials such as the photomultiplier 
tubes and light guides are not shown. 

Mylar. The CS was a 2.5 cm diameter, 0.2 cm thick 
(mass 4.6g) pressed powder disc of Gd 2 Ti 2 07, a geo- 
metrically frustrated antiferromagnetic pyrochlore. At 
room temperature Gd2Ti207 is a paramagnet in which 
fluctuating magnetic moments are known to depolarize 
muon spins within a few microseconds, with an expo- 
nential form of P M (£)i2£ Note that no grease or glue was 
used in the assembly of the targets. Sets A-D used tar- 
gets of pure CS or pure metal, before we realized the im- 
portance of maintaining constant material thickness in 
the decay positron path when making systematic com- 
parisons of decay asymmetry. We anticipated that de- 
termining the fraction of muons stopping in the trigger 
scintillator (TM2) would be a dominant source of sys- 
tematic uncertainty, thus for sets I-K the thickness of 
the trigger scintillator was altered. 



V. ANALYSIS 

When analyzed with a Fourier transform algorithm, 
the data showed time structures with periods of 43.37 ns 
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TABLE I. Data sets for the /i + SR experiment. Sets C and D 
had a low muon rate. Set D had the DC separator on a high 
setting. In the target column, the upstream material of the 
reversible target is listed first. 

Data set Target Scintillator Duration /i + counts 

thickness (hours) Ag mask sample 
(/xm) (xlO 9 ) (xlO 9 ) 



A 
B 
C 
D 
E 
F 
G 
H 
I 
J 
K 



CS^ 
Al 
Al 
Al 

CS+A1 
Al+CS 
Ag+CS 
CS+Ag 
CS+Ag 
CS+Ag 
Ag+CS 



254 
254 
254 
254 
254 
254 
254 
254 
508 
127 
127 



14.1 
8.0 
15.7 
5.7 
24.0 
40.1 
32.4 
22.8 
18.0 
22.2 
13.3 



1.2 
0.7 
0.7 
0.3 
1.7 
2.9 
2.5 
1.8 
1.3 
1.5 
1.0 



0.8 
0.5 
0.5 
0.2 
1.2 
2.0 
1.7 
1.2 
0.9 
1.0 
0.7 



CS: Gd2Ti207 calibration sample 



and 3.68 ns, corresponding to the period of the TRIUMF 
cyclotron and the precession frequency of the muon trans- 
verse polarization components in the 2.0 T magnetic field. 
To minimize the impact of possible effects, we rebinned 
our data using time bins of width 89.84 ns starting 10 ns 
after the muon stop time, and we later verified that our 
results are robust to the binning used. 

We simultaneously fit the B and F decay count his- 
tograms with standard coupled /i + SR equations, 



n b (i) = b b + N e-V T » [1 + A b P^ (t)] 
n f (t) = b f +rN e- t ^[l + A f P^t)} 



(5) 
(6) 



where n b and n f are the number of B and F counts, b b and 
bf are the backgrounds in each counter, r M = 2.19703 /is 
is the muon lifetime, Nq is the histogram normalization, 
A b and Af are the empirical asymmetries that depend 
only on the counter geometry and the material through 
which the decay positrons pass. The parameter r is in- 
troduced to account for differences between the B and F 
counters' energy and angular acceptance. The CS sample 
was fit first, using three forms for the depolarization: a 
single exponential, P M (t) = P/j,(0) exp(— At), a sum of two 
exponentials, and a power law, P/j,(t) — exp(— at p ). The 
data were consistent with a single exponential relaxation, 
so this form is assumed in the following. (See Fig. [2j) 

The backgrounds in each counter were determined us- 
ing the 1 /zs of pre-trigger data, and also by leaving the 
background as a free fit parameter. We found that the 
two measurements differed by more than three standard 
devations, probably due to contamination of the back- 
ground from long-lived muons. Since the backgrounds 
could be well determined using the decay data and had a 
weak correlation with the other fit parameters, we chose 
to fit the background level with a free parameter. 
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FIG. 2. Depolarization for the CS (Gd 2 Ti 2 7 ). (a) rela- 
tive polarization and (b) normalized residuals. The fit is 
P M (t)/P M (0) = exp(— \t) to data from set E, with A = 
(0.860 ± 0.004) /xs" 1 and the x 2 / ndf = 285.5/302. 



VI. RESULTS AND SYSTEMATIC 
UNCERTAINTIES 

The parameters from fits to data for the Ag mask, 
which only required a trigger in TM1, were affected by 
changes in the target and running conditions. This de- 
pendence was seen dominantly in the values of the pa- 
rameter r. The value of A was somewhat sensitive to 
the target, but the average measured value of A mas k = 
(1.07 ± 0.07)ms _1 agrees with the result for the high pu- 
rity Ag given below. Systematic errors were not evalu- 
ated for the Ag mask. 

We used the data from sets I and J to determine the 
relaxation rate for muons stopping in the trigger scintilla- 
tor. Initially the data with a 127 /mi thickness scintillator 
were fit assuming that all muons passed through the scin- 
tillator and stopped in the CS. With the CS relaxation 
fixed, we fit a sum of exponentials to the data with scin- 
tillator thickness 508 /im, which provided an estimate of 
the relaxation within the scintillator itself. This process 
was iterated, and we found that a single exponential is 
adequate to describe the scintillator depolarization, with 
86.0+0.3% of the muons stopping in the 508 /im thickness 
scintillator and A scint = (0.0132 + 0.0008) fis^ 1 . Then the 
CS relaxation is A C s = (0.866 ± 0.008) /is" 1 . With the 
relaxation rate for the scintillator component fixed, an 
attempt was made to determine the fraction of muons 
stopping in the other two scintillators. For the 127 /jm 
thickness scintillator we find (4 + 9)%, and for the 254 /jm 
scintillator we find (0 ± 10)% for set E and (8 ± 8)% for 
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set H. These results are clearly imprecise, so we resorted 
to an evaluation of the fraction using a SRIM 34 simula- 
tion of the muon stopping position distribution. Based 
on estimates of the M20B beam line mean momentum 
(between 29.2 and 28.7 MeV/c) and momentum resolu- 
tion (between 1 and 2%), and of materials in the path of 
the stopping muons, the results of the simulation showed 
that (3 ± 3)% of the muons could have stopped in the 
254 /im scintillator, and no more than 0.2% in the 127 
scintillator. We correct our result assuming that 3% of 
muons stop and depolarize in the 254 [im TM2 trigger 
scintillator rather than the metal target, and assign a 
large systematic uncertainty that allows for between 1% 
and 6% of muons stopping in the scintillator. 

The empirical coefficients (Af, and Af) from the CS 
analysis were determined with statistical precision < 
0.3% and were constrained to be Af, = 0.185 and Af = 
—0.238 while fitting the asymmetry data for the metals. 
There is a negligible contribution to the systematic un- 
certainty from the fixed Af, and Af values; a change of 
10% in either of these parameters altered the relaxation 
rate by at most 0.12 ins -1 . The r parameter was left 
free since it was found to be highly sensitive to the exact 
placement of the target. This is shown in Fig. |3l which 
also includes the results for the relaxation parameter A 
in Eq. gj). 
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the background. Also, the sets with different muon rates 
and DC separator settings showed no evidence of an in- 
consistent relaxation rate. Therefore no systematic un- 
certainties are assigned for these effects. 

The uncertainty in the fraction of muons stopping in 
the 254 /im scintillator is much larger than anticipated, 
and dominates all others. Unfortunately we cannot re- 
duce this uncertainty, since a precise measurement of the 
M20B mean momentum and resolution was not taken 
with our bcamline settings. 

Our results for the 254 /jm scintillator are then 

Aaj = 1.3 ± 0.2 (stat.) ± 0.3 (syst.) ms _1 and (7) 
A Ag = 0.9 ± 0.2 (stat.) ± 0.2 (syst.) ms _1 . (8) 

A single exponential fit is shown for an aluminum run 
in Fig. [4] and for a silver run in Fig. [5] There is no 
evidence that anything beyond a single exponential is 
needed to describe the depolarization. There is no ob- 
servable fast depolarization component below 1 /is. 
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FIG. 4. Aluminum relative polarization vs time, (a) full time 
range, (b) first 6 /xs, and (c) normalized residuals from fit. 
The fit is P M (f)/P M (0) = exp(-At) to one run of set F, with 
A = (1.5 ± 0.5) ms" 1 and the xV ndf = 280.6/304. 



FIG. 3. Results for fits to the metal samples, (a) A for Al, 
(b) A for Ag, (c) r for Al, and (d) r for Ag. Sets F and G had 
multiple runs. The parameter r is highly sensitive to target 
placement, which occurred at the beginning of each set. 

A significant time structure of the background would 
give rise to a systematic uncertainty. Because the confi- 
dence levels for the fits were reasonable, except for one 
set that still produced a consistent relaxation rate, there 
was no need to investigate a possible time structure of 



For the single run with the thin (127 /im) scintillator, 
the contribution from scintillator stops is negligible, but 
the statistical uncertainty dominates so that 

A Ag = 1.2 ±0.4 (stat.) ms" 1 . (9) 

Combining the two results for silver yields 

A Ag = 1.0 ± 0.2 (stat.) ± 0.2 (syst.) ms -1 . (10) 
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FIG. 5. Silver relative polarization vs time, (a) full time 
range, (b) first 6 fis, and (c) normalized residuals from fit. 
The fit is Pfi(t) / -P M (0) = exp(— Xt) to one run of set G, with 
A = (1.5 ± 0.5) ms" 1 and the x7 ndf = 300.2/304. 



VII. CONCLUSIONS 

The measured relaxation rates for silver and aluminum 
differ by less than a factor of two, yet the nuclear dipole 
moments differ by a factor of 35, providing further evi- 
dence that the depolarization is not from nuclear dipole 
moments. These results are consistent with those from 
the final phase of the TWIST experiment,— but with 
somewhat larger uncertainties. The unique conclusion 
from this experiment is that no additional depolarization 
components exist in the time range 0.010 < t < 1.000 /.is; 
this is very important for the interpretation of TWIST 
data. We are unaware of any experimental technique that 
would allow us to readily study the depolarization below 
10 ns, but we are also unaware of any credible models for 
muon depolarization within the first 10 ns in nonmagnetic 
metals. 
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